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We have studied the spin structure of circular four-electron quantum rings using tunable confine- 
ment potentials. The calculations were done using the exact diagonalization method. Our results 
indicate that ringlike systems can have oscillatory fiips between ferromagnetic and antiferromagnetic 
behaviour as a function of the magnetic field. Furthermore, at constant external magnetic fields 
there were seen similar oscillatory changes between ferromagnetism and antiferromagnetism when 
the system parameters were changed. According to our results, the magnetism of quantum rings 
could be tuned by system parameters. 
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I. INTRODUCTION 

During the last two decades, there has been seen an 
increasing scientific and technological interest in spin re- 
lated phenomena and in possible application of these in 
future data processing, communication and storage.^ 
Quantum dots have been proposed as components in 
few-electron spintronics devices, such as spin filters or 
spin memories,—!^ and in spin-based quantum computa- 
tion devicesj^ 

Existing fabrication techniques allow construction of 
semiconductor quantum rings of nanometer dimensions 
containing only a few electronsi^ii Nanoscopic quantum 
rings are sufficiently small systems to show quantum ef- 
fects, and large enough to be able to trap magnetic flux 
in their interior, when subjected to an experimentally 
reachable magnetic field. The trapping of magnetic flux 
quanta gives rise to interesting effects, such as persis- 
tent currents and other periodic properties related to the 
Aharonov-Bohm effect -^:^ 

The main focus of this paper is on the magnetic proper- 
ties of the quantum rings. The previous study by Kosk- 
inen et al.^ for one type of quantum rings has shown 
that a model of localized charges and antiferromagnetic 
coupling of the nearest-neighbor spins, corresponding 
to an antiferromagnetic Heisenberg model, captures the 
physics of systems they study. This means, for example, 
that when the angular momentum of the system is in- 
creased by making the magnetic field stronger, the change 
of the total spin of the system follows values found using 
the antiferromagnetic Heisenberg model. On the other 
hand, a study of a hard- wall quantum dot^ where strong 
electron-electron interaction forces the electron density 
to be ring-shaped, has found a ferromagnetic behavior of 
the system as a function of the magnetic field. Motivated 
by this discrepancy, we study tunable quantum rings and 
aim to identify the underlying physics that leads to dif- 
ferent magnetic properties for the various quantum rings. 

Our results show interestingly that the same quantum 
ring system can show oscillation between ferromagnetic 
and antiferromagnetic behavior as a function of the mag- 
netic field. A similar control of the magnetism can be 



obtained by changing the width or the radius of the quan- 
tum ring. 

II. MODEL AND COMPUTATIONAL METHOD 

We use an effective-mass approximation and model the 
semiconductor quantum rings as two-dimensional sys- 
tems with the Hamiltonian 

where N is the number of electrons, A is the vector po- 
tential of the perpendicular magnetic field B = Buz, V 
is the external confinement potential, m* is the effec- 
tive electron mass, C is the Coulomb constant (normally 
1), and e is the dielectric constant. The electrons are re- 
stricted to the Cartesian xy plane, and the magnetic field 
is normal to the plane. In our calculations, we use the ef- 
fective parameter values of GaAs, namely, to* — O.OGTme 
and e = 12.7. To enhance the spin effects, we have left 
out the Zeeman potential. This can be justified since ex- 
perimentally the gyromagnetic factor, and thereby also 
the Zeeman term, can be made vanishingly small>^ 

To control the spin effects, we use confinement poten- 
tials of the form 

V{r) = ^m*ujy + exp(~rVf7'), (2) 

where determines the confinement strength and Vq 
determines the strength of the Gaussian perturbation. 
We set hjjjQ = 5 meV and cr = 2 a^, where ajj « 10.03 
nm is the effective Bohr radius of GaAs. By tuning the 
Gaussian perturbation stronger, we can make the system 
more and more ring-shaped. In addition, we have used a 
parabolic ring potential given by 

V{r)^\m*ojl{r~r,f (3) 

to model a narrow quantum ring. 
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In the exact diagonalization (ED) calculations, our 
single-particle states are the one-body eigenstates of the 
Hamiltonian H. In cases when the confinement poten- 
tial V is parabolic, the single-particle eigenstates are 
the Fock-Darwin states. For systems with perturbed 
parabolic or parabolic ring potentials, we calculate the 
single-electron states as linear combinations of Fock- 
Darwin states, all with the same angular momentum, 
corresponding to the 15 lowest energy values. The ED 
calculations take into account Landau-level mixing, as 
the 25 lowest-lying eigenstates, irrespective of Landau 
level, are included in the computational procedure. The 
ground state eigenvalues of the Hamiltonian matrix of 
the ED method were obtained using the Lanczos diago- 
nalization method, and the interaction matrix elements 
were calculated using numerical integration. The princi- 
ples of the ED method are described in Ref. [Tsl . 

When the electrons in a quantum ring become suf- 
ficiently localized, charge and spin excitations sepa- 
rate from each otherJ^ This phenomenon is a non- 
perturbative effect due to the strong correlations in quan- 
tum ring systems at high magnetic fields. The qualita- 
tive behavior of the many-particle spectrum of a quasi- 
one-dimcnsional system has been described by the lattice 
Hamiltonian^ii^ii^ 



TABLE I: Total angular momenta L and their correspond- 
ing spin values 5* for the ferromagnetic and antiferromagnetic 
ground states of the Heisenberg model of a four-electron ring. 
Spin states corresponding to higher angular momenta are ob- 
tained by utilizing the periodicity of the spin sequence. 
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iJ = J^S, •S, + — L2 + ^?ia;„na, (4) 
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where the first term is a Heisenberg Hamiltonian that 
models the spin degrees of freedom, the second term is a 
rigid rotation of the system, and the last term includes 
vibrational modes of the localized electrons. The pa- 
rameters of the model are the nearest-neighbor coupling 
constant between the spins J, the total moment of inertia 
/, and the vibration frequency LOa ■ The angular momen- 
tum of the ring is given by L and Ua is the number of 
excitation quanta of the vibrational mode. In an antifer- 
romagnetic system the coupling is such that J > 0, and 
in a ferromagnetic system J < 0. 

The total angular momentum-spin pairs of the low- 
est eigenstates for the effective Hamiltonian Q can be 
calculated using exact diagonalization, or using group- 
theoretical methods.^ In TablelH we have given the angu- 
lar momentum and spin states of the ferromagnetic and 
antiferromagnetic systems of four particles, as calculated 
by the exact diagonalization. The spin sequence is TV- 
periodic as a function of the angular momentum, where 
N is the number of electrons. One can see that for four 
particles, the antiferromagnetic behavior is manifested 
by the three consecutive ground states with spin equal 
to one, and one signature of the ferromagnetic model is 
the occurrence of ground state with spin equal to two. 
One should also note that there is one ground state with 
total spin equal to zero for both versions of the Heisen- 
berg model. 



FIG. 1: (Color online) Ground state spin polarization 
as a function of the magnetic field B for a four-electron 
parabolic quantum dot with Gaussian perturbations of dif- 
ferent strengths Vq. The spin states are (red/middle grey), 
1 (green/light grey), and 2 (dark grey). The upper integers 
denote the number of central vortices found in each of the dif- 
ferent states, and the integers surrounded by a ring denote the 
corresponding angular momenta. All states at zero magnetic 
field have angular momentum 0. If two subsequent states are 
separated by a solid line, the angular momentum increases by 
one when moving to the state at right, and if they are sepa- 
rated by a dashed line, the angular momentum increases by 
two. 



III. RESULTS 
A. Gaussian perturbation ring 

We start with the results for the parabolic dot that has 
been perturbed with the Gaussian potential in the center 
to form a quantum ring. In Fig. [U we have plotted spin 
ground states as a function of the magnetic field B for 
four-electron perturbed parabolic dots with perturbation 
strengths Vo = 0, 1, 2, and 3 Ha* (1 Ha* w 11.3036 meV). 
For all these perturbation strengths, the radius of the 
quantum ring is around 2 a*g. For Vq = 0, corresponding 
to a pure parabolic quantum dot, the states with L — 6 
and 18 are the only fully spin-polarized states (dark gray 
regions in phase diagram), as also found in Rcf. 14 in 
the lowest Landau level approximation. The state with 
L = 6 is called the maximum-density droplet (MDD),— 
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corresponding to the fractional quantum Hall states with 
filling fraction ly — 1, and the state with L — 18 corre- 
sponds to ly ~ 1/3.— Between these two states the total 
spin has values of zero and one, but these do not fol- 
low the predictions from the antiferromagnetic Hcisen- 
berg model, most clearly seen from the fact that there 
are no regions with three consecutive ground states with 
spin equal to one. This is not surprising, as we expect the 
Heisenberg model to be relevant only in the limit where 
electrons are strongly in a ring-shaped confinement. 

Next we turn on the Gaussian perturbation, and one 
can see that the spin-polarized state with i = 6 at Vq = 
splits into several spin-polarized states with different an- 
gular momentum. In addition, between these S ~ 2 
states, the total spins follow the ferromagnetic Heisen- 
berg spins given in Table [H The only exception on this 
rule is found at Vb = 2 around _B = 11 T, where S* = 1 is 
found instead 0. 

Outside these ferromagnetic regions at stronger mag- 
netic fields, one can see antiferromagnetic behavior. The 
antiferromagnetic phase can be identified by the three 
S — 1 states in a row, separated by the 5 = state. 
In addition, antiferromagnetic behavior could develop at 
weak magnetic fields, as at Vb = 3 there are already three 
consecutive states with 5 = 1. As will be shown later 
in this paper, when the ring is sufficiently narrow there 
is indeed found an antiferromagnetic region at magnetic 
fields lower than the first ferromagnetic region. 

In a previous studyiS, of a hard-wall quantum dot with 
a radius R — 5 a*g and maximum of the electron density 
around 3 a|j, the spin structure was ferromagnetic for 
angular momenta L between 6 and 18. Such a magnetic 
structure is almost the same as that found for the per- 
turbed parabolic system with perturbations Vb = 1 and 
2. The magnetic structure of the wide quantum ring is 
thus for angular momenta 6 < i < 18 similar to that of 
a hard- wall quantum dot. 



B. Vortex structure 

From the phase diagram shown in Fig. [U one can also 
see that the state with angular momentum L = 6 changes 
spin when the Gaussian impurity is made stronger. This 
change and many other details of the phase diagram can 
be understood by studying the vortex structure of the 
many-body wave function. 

When subjected to an external magnetic field, the elec- 
trons in a quantum dot are forced to rotate. If the ro- 
tation is sufficiently strong, vortices are formed in the 
electron liquidji^ One can consider the vortices as quasi- 
particles that are holes in the occupied Fermi see.:^ The 
vortices are seen as zero points in the conditional electron 
density, where the phase of the wave function changes as 
a multiple of 27r for a path enclosing the node. As shown 
in the fractional quantum Hall effect (FQHE) theory,^° 
the formation of vortices is a result of quantization of 
the magnetic field, and each vortex is associated with an 
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FIG. 2: Conditional wave functions consisting of the charge 
density, plotted with contour lines (logarithmic scale), and 
the phase of the wave function, illustrated by the gray-scale 
shading. The fixed spin up electrons are marked with black 
triangles pointing up, and the spin down electrons with white 
triangles pointing down. The probe particle is marked with 
a larger triangle. The clusters consisting of central vortices 
(nodes), which can be considered as giant vortices, are de- 
noted by circular arrows. 

integer number of Dirac flux quanta. 

To analyze the nodal structure of the perturbed 
quantum dot, we have calculated conditional wave 
functions defined for an iV-particle system as 

where r is the position of the moved particle 1 , and r* is 
the most probable position of particle i that are found by 
maximizing the total electron density The phase 9 

is obtained from the relation V'c(r) = |'!/'c(r)| exp(i6'(r)) 
and the conditional electron density is defined as |'(/'c(r)p. 

Now returning to the question related to the spin of the 
L = 6 state, a starting point for this analysis is the con- 



4 



ditional wave function of the parabolic dot shown in Fig. 

In this, one can see a Pauh vortex on top of each 
electron, except the probe particle on right. The vortices 
are shown by the discontinuous jumps in the gray-scale 
from white to black when the electrons arc circulated in 
a clockwise fashion. The fact that the vortex number is 
the same as the electron number shows that the state is 
a finite-size example of a quantum Hall state with filling 
fraction v ^ \. 

For Vb = 1, shown in Fig. [Hg), the spin is still the 
same and the conditional wave function is nearly identi- 
cal to the Vb = state. However, for Vq — 2 the spin has 
changed from 2 to 0, and the conditional wave function 
in Fig. [2{d) looks completely different. There is still one 
Pauli vortex on top of the left-most electron, but the two 
other electrons have switched spins. In addition, there is 
now one vortex that is located at the center of the sys- 
tem. It turns out that one needs to combine two Pauli 
vortices in order to make one vortex at the center of the 
dot. From this data, one can understand the change in 
the total spin as follows: The system places the vortices 
so that it minimizes the total energy. Without the Gaus- 
sian impurity at the center of the dot, it is energetically 
favorable to place the vortices on top of the electrons to 
reduce the Coulomb repulsion. When the potential at 
the center of the dot is raised by the Gaussian impurity, 
at some strength it is more favorable to reduce the prob- 
ability of the electrons to be at the center of the dot by 
placing a vortex there. 

One can do a similar analysis for the L — 1% state, 
which for Vb = corresponds to the ;/ = 1/3 state. The 
correspondence can be seen by the three vortices bound 
to each electron, as shown in Fig. [H^k). The reason for 
the small separation of the vortices is due to long-range 
nature of the interaction. For Vb = 1 shown in Fig. [2Ih), 
there are already two vortices at the center of the dot, and 
to create these, four Pauli vortices have vanished from the 
system. Due to this, the spin has again changed, and the 
opposite spins have only two vortices bound to them. As 
a general rule, opposite spin electrons can have even vor- 
tices bound to them, and the antisymmetry requirement 
of same spin electrons forces the vortex number to be 
odd. The vortex structure of the 5 = states is similar 
to that oia,v = 2/3 Halperin statej2ii^i22 apart from the 
vortices at the center of the dot. Going to Vb = 2 shown 
in Fig. [2Ie), a third vortex has appeared at the center, 
and now again the total spin of the system has changed. 
Finally, at Vb = 3 shown in Fig. Hfb) the conditional 
wave function is nearly identical to the Vb = 2 case. 

The same trend can be found at even larger values of 
the angular momentum, and as a final example, we show 
the data for the L = 30 states in Figs. [2ljc),(f), and 
(i). The main difference to the previous examples is that 
when the number of central vortices grows, also the area 
they occupy seems to be larger. 

Analyzing the vortex structure of the other states in 
detail enables us to label the central vortex numbers of 
the ground states. These numbers are given in Fig. [1] 



for the most interesting states. One can see that in gen- 
eral, the number of central vortices grows as a function 
of the magnetic field. However, at each fixed Gaussian 
impurity strength, there is a point where it is energeti- 
cally favorable to add Pauli vortices instead of the central 
vortices, and at this point the number of central vortices 
is constant although the angular momentum increases. 
At these points, the ferromagnetic behavior is changed 
to antiferromagnetic. Somewhat similar transitions are 
seen at the magnetic fields below the first ferromagnetic 
states. 

It is also interesting to analyze the ground states at 
magnetic fields around 6-7 T. At Vb = 0, the ground state 
is the 5* = 2 state with L — 6 corresponding to — 1. 
When the Gaussian impurity is made stronger, there are 
more and more vortices at the center of the system, and 
the angular momentum is increased. However, the spin 
of the ground state does not change. This shows that 
the V = 1 state is in some sense stable even for this small 
particle number when one pierces it with three fluxes at 
the center. By this we mean that the vortex structure of 
the conditional wave function is the same apart from the 
central vortices, see Figs.[IJb) and (j). 

C. Narrow quantum ring 

Based on the analysis presented above, one would ex- 
pect to find an antiferromagnetic region at low magnetic 
fields, if the ring is made sufficiently narrow (see the re- 
gion with S = 2-4 T and Vb = 3 Ha* of Fig. 1). To 
investigate this, we switch the confinement potential to 

V{r) = ^m*^l{r-rof, (6) 

where the confinement strength is taken to be 40 meV 
and the radius So^/tt 2.5a*g. This system is a very 
narrow ring with approximately the same radius as the 
ring studied above. 

In Fig. [21 the total angular momentum and spin of 
the ground state is plotted for the narrow ring. We have 
varied the Coulomb constant between the values C = 
0.1, 1, and 10. In reality, the interaction strength can be 
changed by changing the quantum ring radius. ^2, The case 
with C = 1, shown in Fig. OJb), corresponds to the natu- 
ral interaction strength, and serves as basis for the com- 
parison with the ring used above. For this case, the total 
spin has values zero and one, and above _B = 4 T there 
are three states with spin one between each ground state 
with spin equal to zero, in agreement with the antiferro- 
magnetic Heisenberg model. Thus, as expected there is 
an antiferromagnetic region at low magnetic fields when 
the ring is sufficiently narrow. 

When the Coulomb interaction is made weak by setting 
C = 0.1, two of the three consecutive S = 1 states are 
no longer ground states as shown in Fig. El^a). Now, 
although the spin has values and 1, the system does 
not behave as the antiferromagnetic model. This shows 



5 




B(T) 



B(T) 



B(T) 



10 12 14 



FIG. 3: (Color online) Total angular momentum L and spin 5* of the ground states given as functions of the magnetic field B 
for a four-electron parabolic ring with 40 meV confinement strength and radius %a*g/iT, where a*g is the effective Bohr radius. 
The Coulomb constant is (a) C = 0.1, (b) C = 1, and (c) C = 10. The crosses denote states with magic angular momenta 
Lk ~ iMDD + kN , where k G {0, 1, 2, . . . }, Lmdd — {N — l)N/2 is the angular momentum of the maximum-density droplet 
(MDD), and A'' is the number of electrons. The state with angular momentum Lmdd is marked with a cross surrounded by a 
circle. 



that the localization induced by the strong interaction is 
a necessity for antiferromagnetism. 

On the other hand, when we make the interaction 
stronger by setting C — 10, one can see that fully spin 
polarized ground states are found. In addition, the mag- 
netic field range of each ground state is nearly the same. 
In the magnetic field region between 6 and 10 T, the 
ferromagnetism is not as pure as the antiferromagnetism 
for the C = 1 case, as now the ground states with S = 
are missing. Also, the system is antiferromagnetic upto 
B K 6 T. We conclude that a transition between anti- 
ferromagnetic and ferromagnetic behaviour also can be 
controlled by the Coulomb interaction strength, or phys- 
ically by tuning the radius of the quantum ring. 



the magnetic phase is determined by the magnetic field 
region, the ring width and the radius of the system. We 
found that the magnetic phase changes oscillatory as a 
function of the different system parameters. Knowledge 
of the nontrivial magnetic phase structure could be very 
fruitful for experiments, as the observed phenomena 
opens the possibility for tunability of the magnetism 
by changing system parameters. For example, in an 
experiment, a metallic electrode with adjustable voltage 
could be used to achieve the control we had in our 
calculations when using the Gaussian impurity. 
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